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Determinations of the gluon propagator in the continuum and in lattice simulations 
are compared. A systematic truncation procedure for the quark Dyson-Schwinger 
and bound state Bethe-Salpeter equations is described. The procedure ensures 
the flavour-octet axial-vector Ward identity is satisfied order-by-order, thereby 
guaranteeing the preservation of Goldstone's theorem; and identifies a mechanism 
that simultaneously ensures the absence of diquarks in QCD and their presence in 
QCD^':"^, where the colour singlet diquark is the "baryon" of the theory. 



The Dyson-Schwinger Equations (DSEs) provide apPoincare covariant, 
nonperturbative, continuum approach to studying QCD^I in which the fun- 
damental quantities arc the Schwingcr functions (Euchdean Green functions). 
The weak-couphng expansion of the DSE for a given Schwinger function gen- 
erates its perturbative series. However, it is the nonperturbative nature of the 
DSEs that is most interesting because it entails that they provide an ideal 
framework for the study of confinement, dynamical chiral symmetry breaking 
(DCSB) and the identification of observable effects of bound state substructure 
in interactions involving hadrons. Recent applications of the DSEs have in- 
cluded the study of|L meson spectroscopy;^ tt-tt scattering^ Lu-p mixing and the 
uj-p mass splitting^ the electromagnetic fomi, factors of charged ancLneutral 
pions and kaonsp the anomalous j*Tr^ — > 7 au, 777* — > tlzeJj and 77rpLl transi- 
tion form factors; the electroproduction of vector mesonsjij and deconfinement 
and chiral symmetry restoration in finite temperature QCD.Eil 

In present phenomenological applications the two most used DSEs are 
the gap equation, which yields the dressed-quark propagator, and the Bethe- 
Salpeter equation (BSE) for two-body bound states, which yields the mass 
and bound state amplitude. The primary element of the kernels of these 
equations is the two-point gluon Schwinger function (dressed-gluon propaga- 
tor), D^^{k) = {5f^^ — kf^ky/k"^) Dxik"^)- The two-loop perturbative result for 
Dxik"^) is quantitatively rehable for fc^ > 1 — 2 GeV^, however, for k"^ < 1 GeV^ 
nonperturbative methods are required to calculate D-rik^)- 

Nonperturbative studies of the DSE for the gluon vacuum polarisation 
indicate a strong enhancement of Dxik^) for fc^ < 1 GeV^, with qualitative 
agreement that Dxik^) exhibits a regularised ipficared (IR) singularity, which is 
often characterised as a regularisation of 1/A:^J13'EJ This is illustrated in Fig. 1. 
The regularisation is crucial since Dxik^) appears in integrands sampled in 
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domains containing k'^ = 0. This behaviour is consistent with confinement 
because D^y(k) thus described: 1) doesn't have a Lehmann representation and 
therefore no asymptotic gluon excitation ia_associated with it; and 2) provides 
for area-law behaviour of the Wilson loopO 
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Figure 1: A comparison of possible behavioufs-iof the gluon propagator in the IR. A « 
0.44 GeV in the IR-enhanced gluon propagator 113 b = 0.34 GeV in the IR-vanishmg form:M 
M = 0.13 GeV in the massive-boson form.Ej The error-bars on the lattice pointsEj are only 
my crude estimate of the statistical error associated with 25 configurations. 

The possibility that Drik^ — 0) might be finite or zero has been explored 
and found to be phenomenologically untenable because such a form does not 
provide for quark confinement and oiih^ leads to dynamical chiral symmetry 
breaking under special circumstances .Ej Furthermom-iecent analyses indicate 
that the gluon DSE does not admit such solutions.E3E3 

One can compare these observations with numfiiical simulations in_which 
the Landau gauge gluon propagator is extracted. E3£3 In Fig. 1 resultal3 of a 
simulation with 25 gauge-configurations at /3 = 6.0 on a 16'^ x 40-lattice are 
plotted. It is reported that a x^"fit to all the collected data favours an IR- 
vanishing form of Dxik^), although a massive boson propagator could not be 
eliminated from consideration. Consistent with this uncertainty, one notes in 
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Fig. 1 that it is only the first two data points that begin tojdistinguish between 
the forms plotted. This highlights a recognised difficultytj in extracting the 
small-A:^ behaviour of Drik^) in lattice simulations. Finite lattice size entails 
that few data points can be collected at small-fc^ . 

The true magnitude of this difficulty is greater, as illustrated in a study 
of 500 configurations at /3 = 6.0 on a 24^^ x 40 lattice,E£l which finds that 
Dxik^ = 0) is finite and nonzero. It is argued therein that their first 3 data 
points, at k'^ < 0.34 GeV^, should not be included in the fit: without them 
w 1 and it is stable; including them dramatically degrades the fit-quality. 
(A similar effect is observed if data at large-fc^ is included, which are affected 
by finite-spacing artifacts. It is argued that only the data in a window at 
intermediate fc^ should be included in the fitting procedure.) In the context 
of Fig. 1, such considerations suggest that at least the first two lattice data 
points can be neglected, in which case no curve is favoured over another. 

In lattice simulations, the finite volume provides an intrinsic IE_ cutoff, 
which necessarily entails that the gluon propagator is finite at — 0X3 There 
is a mass-scale, M, associated with this finite value. A first estimate shows that 
M remains finite and nonzero in the infinite volume limit. However, the infinite 
volurue limit alone does not indicate the behaviour of M in the continuum 
limit ,E3 which corresponds to ^ — > oo and /3 — > oo. There is presently no data 
that can provide an indication of the value of M in the continuuna limit. One 
notes that should M in the continuum limit then this studyli£l would yield 
an IR-enhanced gluon propagator. 

For the present the most reliable determination of the qualitative behaviour 
of Drik^) in the IR is provided by the DSE studies. One result of the IR- 
enhanced form of Dxik^) is DCSB without fine tuning. The quark propagator 
is commonly written S{p)~^ = 17 ■ pAlp^) + B{p^) and the quark-DSE is 

tj-p[A{p')-l]+B{p')=m+^ I ^g^D,,{p-k)j,S{k)n{k,p), (1) 

where T^(k,p) is the dressed-quark-gluon vertex. At any finite order in pertur- 
bation theory the dressed-quark mass function, B{p^), vanishes if the renor- 
malised quark mass is zero; i.e., in the chiral limit. The quark condensate is 
nonzero if-and-only-if B[p^) is nonzero. With an IR-enhanced gluon propaga- 
tor (0) necessarily admits BijP') nonzero, even in the chiral limit; hence one 
has DCSB. This result does not depend on the exact form of the Ansatz used 
to describe the quark-gluon vertexE Another result is quark confinement. S{p) 
obtained as a solution of (|^) with an IR-enhanced gluon propagator is itself en- 
hanced in the vicinity of A:^ = and does not have a Lehmann representation. 
This entails that no asymptotic quark excitation is present in the spectrum. 
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All studies of meson and quark-quark (diquark) bound states to date have 
used the rainbow-ladder truncation of the quark-DSE and meson/diquark- 
BSE. The rainbow quark-DSE has rg(fc,p) = 7^ in (|l|) and the meson ladder- 
BSE is 

Tip;P) = -J ^g'D,4p^q)j,^Siq+)riq;P)S{q^h,,^, (2) 

where q± = q i P/2, p is the relative quark- antiquark momentum and P 
is the total momentum of the meson. This pairing has been phenomeno- 
logically successful for ground state flavour-octet [(8)/] pseudoscalar, vector 
and axial- vectoiLJnesons, primarily because it is a Goldstone theorem preserv- 
ing truncation.lij This can be seen heuristically by substituting T]^ {p; P) = 
T'j5[iE^ip; P)+l- PF^ip; P)] into (|) to obtain 

with a similar but more complicated equation for F^^. In the chiral limit, 
rainbow-truncation of (|l|) 

Hence, for = 0, the solution of (|) is E{p;P) oc B{p'^), with F^{p;P) ^ 
and completely determined by A{p'^), B{p'^). In this truncation therefore DCSB 
necessarily entails, without fine-tuning, a massless, pseudoscalar bound state 
of a dressed-quark and -antiquark whose bound state amplitude is completely 
determined by the nonperturbative, dressed-quark propagator. (This result 
persists if the remaining Dirac amplitudes are retained in F^.) 

Underpinning this result is the fact that the rainbow-ladder truncation is 
a (8) / axial- vector Ward identity preserving truncation. The inhomogeneous 
ladder-BSE for the (8) / axial-vector vertex is 

lT%p- P) = n57p -\j D^,{p - q) 7p S{q+)lTl{q- P)S{q^h, . (5) 

Contracting both sides with P^ one finds that the chiral limit axial- vector Ward 
identity: —iP^ ^^^{k; P) — S^^{k+) 75 -I- 75 S~^{k-) , is satisfied if-and-only-if 
S(j)) is the solution of the rainbow quark-DSE. (In considering renormalisation 
this heuristic outline acquires some subtleties but the result is qualitatively 
unchanged.) The conclusion is that Goldstone's theorem is manifest in any 
DSE truncation scheme that preserves the (8) / axial- vector Ward identity. 
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One systematic procedure for constructing such a scheme is based on a 
skeleton graph expansion of the dressed quark-gluon vertex in (|i|)c!j In this 
skeleton expansion every line and vertex is considered to be fully dressed 
except the quark-gluon vertex, which is bare. It is easiest explained via il- 
lustration. The first term, 0{g^), yields the rainbow quark-DSE. Consider 
the integrand in (Q) with rf,{k,p) = 7^: the replacement TZ = ^^S{k)^i, — > 
^^S{k+)T{k\ P)S{k^)ji^ yields the ladder kernel for the meson BSE. This pair 
of equations is a (8) / axial-vector Ward identity preserving truncation. 

In this expansion the O(g^) contribution to the quark-DSE is 

+ iVai3j{kJ,p)Df,a{p - k) D,^p{k - I) Dpj{l - p)-ff,S{k)j,^S{l)j^^ , (6) 

where Vap-y{k,l,p) is the dressed 3-gluon vertex. Performing the replacement 
TZ sequentially at the site of each S in (^) yields, from line one, 2 quark-gluon 
vertex correction terms and 1 crossed-box term, and from line two, 2 3-gluon 
vertex terms. These are the 5 contributions to the meson BSE kernel that are 
sufficient and necessary at this order to ensure that the (8) / axial-vector Ward 
identity is satisfied, which ensures that Goldstone's theorem is preserved at this 
order without fine tuning. This procedure can be continued order- by-order. 

The ladder BSE is purely attractive in the colour-singlet [(l)c] pseudoscalar 
meson channel. Repulsive terms only occur at 0(17^): the quark-gluon-vertex 
correction terms obtained from (|^) via TZ are attractive; the crossed box term 
is repulsive; one of the S-gluon-veitex contributions is attractive, the other 
repulsive. A simple heuristic studycj shows that, at O(g^), the attractive con- 
tributions almost completely cancel the repulsive terms. The terms themselves 
are not small but their sum is. A persistence of this cancellation order-by-order 
explains the success of the rainbow-ladder DSE-BSE pairing for the (8) / pseu- 
doscalar, vector and axial- vector mesons. In the scalar sector all 0(5'*) are 
repulsive and there is no cancellation. This explains the failure of the rainbow- 
ladder pairing for scalar mesons: they are not simply ladder bound states of a 
dressed-quark and -antiquark. 

Givep^the meson BSE it is straightforward to obtain the analogous diquark 
equationcj If a solution of this equation exists then the QCD spectrum contains 
coloured bound states, either colour-antitriplet or -sextet; colour confinement 
entails that no such solutions should exist. The diquark ladder-BSE is 



(p; P)^-J^y D,, ip - (?) 7p Y siq+)rD (<z; p) 



(7) 
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where T indicates matrix-transpose. (|^) binds colour-antitriplet diquark bound 
states in QCD.I3 (There are no solutions of (Q) in the colour-sextet sector just 
as there are no colour-octet solutions of This failing is due to the purely 
attractive nature of the kernel described above. 

It is in this connection that the repulsive terms appearing at O(g^) and 
higher are significant. The algebra of 5'C/(3)coiour entails that the two 0(5'*) 
3-gluon vertex diagrams still contribute with opposite signs in the diquark 
equation, however, relative to the vertex correction contributions, the repulsive 
0{g'^) crossed-box term in the diquark equation is five times stronger than in 
the meson equation. The repulsive effect of this term is further amplified by 
the IR enhancement of S{p). These effects together act to ensure there is no 
solution of the 0(17^) diquark BSE. The persistence of this effect at higher 
order would explain the absence of diquark bound states in QCD. 

A qualitative check of this mechanism and systematic truncation procedure 
is found in QCD^"^^^, where the (l)c "baryon" is a diquark. A truncation 
procedure that, order-by-order, ensures the absence of coloured diquark bound 
states in QCD should simultaneously ensure the existence of (l)c diquarks 
(baryons) in QCD^-=2. 

The generators of SU{2) arc {T'/2}i^i 3, where r are the Pauli matrices. 
The colour structure of the meson is described by the 2 x 2-identity matrix 
and that of the diquark by ir^. At a given order the meson BSE equation will 
involve strings of the form 

7^„T'"5(fc„) . . .7^^^,r'^+^5(fc,+i)rM((?; F)5(fc,)7M.r^^- ■ • ■ S{h)^^y^ (8) 

while the associated term in the diquark equation will be 

7M.^'"'5(fc«)... (9) 

Defining r£ = TdC^, where C = 7274, and using: rV^[T^]^ = —t'^t^t'^; and 
[7^]^ = -C^j^C, @ becomes 

7^„r'"5(fc„) . ..j,^^y^+^Sik,+^)T%iq;P)Sik,h,y^ . . . S{k,)^^y^ . 

(10) 

This demonstrates that, in QCD^^~'^, Vm{(1\ P) and r£(g; P) satisfy the same 
equation, order- by-order. The truncation scheme and diquark confinement 
mechanism therefore satisfy the constraint indicated above. In fact, one sees 
that in QCD^'^^^ the spectrum of mesons and baryons (diquarks) is identical 
(neglecting electroweak effects). In particular, the existence of a pseudoscalar 
Goldstone boson entails the existence of a massless scalar baryon (diquark). 
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These results are simply a manifestation of the equivalence of the fundamental 
and conjugate representations of SU{2). 

This work was supported by the US Department of Energy, Nuclear 
Physics Division, under contract number W-31-109-ENG-38. 
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